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The Schrodinger equation, in hyperspherical coordinates, is solved 
in closed form for a system of three particles on a line, interacting via 
pair delta functions. This is for the case of equal masses and poten- 
■»^ , tial strengths. The interactions are replaced by appropriate boundary 

1-^ ' conditions. This leads then to requiring the solution of a free-particle 

^~^\ Schrodinger equation subject to these boundary conditions. A general- 

ized Kontorovich - Lebedev transformation is used to write this solution 
as an integral involving a product of Bessel functions and pseudo-Sturmian 
functions. The coefficient of the product is obtained from a three-term 
recurrence relation, derived from the boundary condition. The contours 
of the Kontorovich-Lebedev representation are fixed by the asymptotic 
conditions. The scattering matrix is then derived from the exact solution 
j^ ' of the recurrence relation. The wavefunctions that are obtained are shown 

\—{ to be equivalent to those derived by McGuire. The method can clearly be 

applied to a larger number of particles and hopefully might be useful for 
unequal masses and potentials. 

Introduction 

Three-body systems and processes are of fundamental interest in physics ^. 
One of these, with which a number of us have been concerned, is the recombi- 
nation of three-particles to a dimer plus a free particle, in a many body system 
forming a Bose-Einstein condensate [2]. The condensate is not the lowest state of 
the system, but a metastable state. The 3-body recombination is the dominant 
mechanism for cooling and lowering the overall energy of the system. 

Experimental and theoretical studies have shown that this recombination 
rate depends mainly on the two-body scattering length a [HI S El E] , as the 



C^ 



collision energy is low and the interaction is weak - owing to large interparticle 
distances, and on the bound state energies. 

This would suggest that zero-range potentials (ZRP), defined in terms of the 
scattering length [Jj, 



lini 



1 d{r^) 



rij) dr 



-l/a. (1) 



can be applied to model the interaction between the particles of the condensate. 
It has been shown by Nielsen and Macek using the hidden crossing technique 
that the ZRP describes properly the recombination transition in a system of 
three ''^He atoms 0. Also, Gasaneo and Macek showed that the ZRP gives a 
quite good representation for the adiabatic potential of the same system [H]. 
A closed form solution for a system of three-particles interacting via a ZRP 
has been recently presented by Gasaneo et al [2]. The fragmentation process 
^He2 + ''He — s- ^He-I- "^He-I- ^He was studied and relatively good agreement 
was found when compared with the hidden crossing calculations. 

In this paper, we seek to apply our techniques to a famous model: 3 particles 
in one dimension, subject to pair delta-function interactions. For this model, 
introduced by McGuire JOI, one can obtain exact solutions for the wave func- 
tions, the scattering matrix and the binding energies, in the case of particles 
of identical masses and equally weighted interactions. As such it has been ex- 
tended to a larger number of particles Ql], using Bethe's Ansatz ^^I; and also 
found to be exceedingly useful when used as a test-bed for the development of 
a number of different methods (pertubative, Faddeev, hyperspherical adiabatic, 
etc.) US]. 

Here, we note that using ZRP and Q, in 3-dimensions, leads to the Thomas 
effect and the collapse of the 3-body ground state |1^. However, in one dimen- 
sion, an equation similar to 1^ - with a not the scattering length - provides 
boundary conditions which correctly characterizes the wave functions and re- 
place the use of the (5-function interactions, and should therefore again give us 
exact results. One of these, though, is that the recombination rate, for this 
model, is exactly zero. 

In section II we propose a solution, written in integral form, for the free 
particle Schrodinger equation, written in hyperspherical coordinates. A linear 
combination of free particle solutions can then be found to satisfy the boundary 
condition that we alluded to earher, and thus provide us with the solution of 
the problem with interaction. The requirement that the wave function satisfy 
the boundary conditions leads us to one of the important results of this paper, 
namely that the weight of the free particle solutions, in the integral form, satisfies 
a recurrence relation similar to that obtained in the refs. jH] and jH]. 

In section III the method is appHed to a particular case in which two of 
the particles are bound. It is shown that the recurrence relation, defining the 
coefficient of the free-particle expansion, can be solved in closed form and, thus, 
the scattering matrix is also obtained in a closed form. This allows us to have 
a detailed test of our method. In this section it is also shown that the wave 



function obtained is equivalent to the McGuire plane wave solution, and that 
our expression for the S matrix is the matrix obtained by McGuire, in the 
particular case discussed in this paper. In section IV the relation between the 
hyperspherical adiabatic approach and the present one is discussed. 

In Appendix A, the pseudo-Sturmian functions are derived. In Appendix B, 
the wave function is written as the symmetric wave plane in cartesian coordi- 
nates. 



Exact Integral Representation 

To begin the study of the three identical-particle system (therefore with equal 
masses), consider the center of mass and Jacobi coordinates, 

r ^ ^[Xi+X2+ X3), 

. ^ /!(^--) (^) 

the Xi give us the locations of the 3 particles along the line, see Fig. 1 . 
Using polar coordinates, the 2 Jacobi variables allow us to define, in turn, a 
hyper radius R and an angle as 

Tj^Rcose $,^Rsme (3) 

where —tt < 9 < tt andO < i? < 00. In terms of these coordinates the 
Schrodinger equation for the 'relative' system can be written as 

x7T7 

H^{R,0) = i—)E'i'{R,0) (4) 

where 

1 d ^ d 1 52 A 1 

RdR 

The function C {9) is defined by 



"'--^Rm^'m^R^wrR''^'^- ^'^ 



c{9)^\cY,no-e,). (6) 

where the coefficient c equals (3/7r\/2)(2m/fi^)(7, g being the strength of the 
interactions. This c is negative for attractive interactions and positive for re- 
pulsive ones. The angles 9j equal (2j + l)7r/6. The lines 9 = 9j divide the 
(p, 9) plane in six regions. In each region the order of particles is fixed, so that 
between 9^ and 9^, xi < X2 < X3, etc. A different permutation of particles is 



Figure 1: One of three sets of Jacobi coordinates for the three particles. 



associated to each region. >From now on, we will choose the units such that 
2m = 1 and h = 1. In each sector, we now seek a free particle solution that 
satisfies the boundary condition that will replace the effect of the potential, i.e. 



lini 



1 5*(i?, 



R-iiiR^i 



de 



1 



(7) 



where 6^ = 6 < 6j j = 0,1,. ..,5. In Eq.Q a — {6/ttc) does not depend on 
j because all the strengths of the interactions and all the masses are equal. 
Writing a solution for the free particle system as the product V'/ree (R, 6) = 
6 iy, 9) R^/'^TZy (KR), where K^ — E, leads to the set of free particle equations 



i?2 



i?-i/2z, (Ki?) 

-1 /d^ 



(9^ 
-dR' 



K' 



e(.,.)V^^'' + iJ^(^'') = '^^ 4 



R-^^^Zy (Ki?) 
1 



(8) 



where Zi, (Ki?) = R^'^TZi, (Ki?) is a Bessel function and v a separation constant. 
If for the Q{h',9) functions we choose the pseudo-Sturmian functions S(v,6), 
defined, for fixed v, as the solutions of 
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p{v)C{e) 



S{y,e) 



s{iy,e), 



(9) 



then the functions V'/ree (6*, i?) = S{i', 9)Z„ (Ki?) are solutions of the Schrodinger 
equation Eq.l^J for values of plv) = R. Note that Eq.jsj may be replaced by 
the relation 



subject to the boundary conditions 
"1 I dS 



Siiy,e)^0 



lim 



{,,,9) 



--, J=0,1,, 
a 



(10) 



(11) 



_p{iy)S{:^,0) de 
In the last equation, we assumed that S{i', 9) is symmetric about each line 



We now propose to write the general wave function of the system as a 
Kontorovich-Lebedev transform, in terms of the base functions just discussed, 
that is, as 



*(i?,6')= I dvA{v)S{v,e)Z^{KR), 



(12) 



provided that its derivative satisfies the boundary conditions, i.e. Eq.Q. The 
contour of integration must be chosen so that the wave function has the correct 
asymptotic behaviour. 



Following the reasoning of Gasaneo et al jH], we will now show that the 
boundary conditions, Eq.Q, can be transformed into a recurrence relation for 
A (zy). First, we substitute En. ((T2ll in Eq.Q, and then interchange the order in 
which the integral and the derivative are taken, to obtain for each j 

dS{v,e) 



lim / dvA{v) 

" -♦^J J<; 

— S{iy,e)Z^{KR) 



j^M^R)- 



89 



0. 



(13) 



Second, we use Ea. (jll|l and the identity (2iy/z)Z^{z) = Z^^i{z) — Z^^i{z) to 
transform the equation to 



lim 



diyA{iy)-[-p{i^)/a]S{i^,e) 

V 

x[Z,+i{KR)-Z,^i{KR)] 

+ J- f diyA{iy)S{iy,e)Z^{KR)\ = 0. 



(14) 



We assumed in the previous equation that K = iK and K >0, because we are 
mainly interested in negative energies. By selecting the appropriate contours 
we can now transform the last equation to 



lim 



dv 



A(^-l) 



1 



1 



I/- 1 



p{v- l)S{v-l,t 



-A(v + 1) -p{v + l)S{v + 1, 61) 



K 



A{v)S(y,e) 



Z^{KR) = 



(15) 



Since the set of Bessel functions forms a complete set of basis functions, 
the function within the square brackets should be zero at the limit. We arrive, 
finally, at the recurrence relation that we are looking for. 



B{v-l)p{v-l)S{v~'l,ej) 



(16) 



2v 



-B{u + l)p{v + l)S{v + 1, e,) = —B{v)S{v, Oj) 

K 



where B{v) = A{v)/i>. In the following section, we will apply this approach to 
a particular case of this three body system and show that we can obtain the 
wave function and the 5-matrix. 



2+1 System 

Consider now the case where two of the particles are bound. The wave function 
^{R, 0) can still be written in terms of the Kontorovich-Lebedev representation. 



Ea. lll2ll . The unnormalized angle pseudo-Sturmian function S{v,0), a six-fold 
symmetric function, is defined by the Eas. ljlO|l and ljll|l . As can be seen in 
Appendix A, the function 5'(z/, 9) may written as 



with j — 0, 1, ..., 5, where p{h') satisfies the relation 



J: 



< 



6' 



z/tan(r/-) = — — 
6 (d/ttc) 



(17) 



(18) 



>From the previous section, we can immediately conclude that A{v) satisfies 
the recurrence relation 



yl(i^ + l)sin (iy + 1)- -A(i^-l)sin (u - 1) - 
6 J L 6 



3K 



A{i') cos 



(19) 



Solution of the recurrence relation 

The recurrence relation, displayed in En. lfTfl)! . can be written as 



J{TT/6)iy 



A{iy + l)e'^''/^^ ~A{v~l)e- 



i(ir/6) 



ITTC 



Ail.) 



-i{7r/6)i' 



+A{i. - l)e'('^/6) + ^A{,y) 



= 0. 



(20) 



An inspection, of the solution of the recurrence relation - Eq.(25) in ref. [H], 
leads us to propose a coefficient in the form of the series 



Ail.) 



^~fi^ 



i(Tr/3)i/ _|_ ^gj(7r/3)i^ 



+5^e-*(^/6)^ + 526*^'^/^)'^ + 53 



(21) 



Substituting this expression in En. (l2?1)l . and equating to zero the coefficients of 
exponentials, with different arguments that depend on v, we obtain the following 
values for the parameters: 



S — tan i/3 cot 

\6 / V6 



S3 — — cot — cot 
6 



ih'^) 



(22) 



sin(i^) = z— . 
Consequently, the solution for the coefficient can be written as 

A {v) = e-'^" U-'''/^> + S e("/3)'' + ^3 ) , 



or 



where 



and 



A{v) =2e-'^''\cos{^v + 5\ +, 



S = e 



2i5 



1 TT , 

a — — cot — \ / cot 

2 6 



(J-./3)cot(^+z/?) 



(23) 

(24) 
(25) 

(26) 
(27) 



In the next section, we demonstrate that S represents the scattering matrix 
and, accordingly, 6 the phase shift. We should stress the remarkable fact that 
the S-matrix appears explicitly in the solution of the recurrence relation. In 
the next subsection it is shown that the expression obtained for S in this work 
is equivalent to the formula for the exact symmetric iS-matrix for the 2 + 1 
process, given in Ref. [T7| . 

Asymptotic wave function 

To be specific we will restrict the following discussion to the case of total neg- 
ative energies, and will write K — A/(7rc)2/36 — k^ > 0, in which — (7rc)^/36 
is the two-body bound energy and k^ is the effective energy. Next, we will 
show that the imaginary axis is the appropriate contour to obtain the correct 
asymptotic behaviour of the wave function. Substituting the coefficients A(i/) 
defined in Ens. ll22ll and lf2ijl . the pseudo-Sturmian functions given in En. lfTTIl 
and the modified Bessel functions K^{KR), into En. lfT2|l . as well as choosing 
the imaginary axis as the contour of integration, we find 



X cos 



dv (cosh [{in/S + P) v] - sinh [(i7r/3 + (3) 1^]) 
+S / dv (cosh [{in/ 3 - P) v] + sinh [{in/S - P)v]) 



-^3 / dv (cosh [(3i'] — sinh [(3iy]) 



X cos 



J3)- 



K,{KR). 



(28) 



Note that in the above expression the exponentials in the coefficients have been 
written in terms of hyperboHc functions. The integral over the odd terms van- 
ishes, leaving only the even terms in the integrand. After a trigonometric iden- 
tity, this leads to 



f = 
1 



dv cosh 



dv cosh 



P + ^[e-[J-l]-] v]K,{KR) 
-l3 + i{e-[j + l]'^]\u]K,{KK) 



I dv cosh ( ^P + i{e-[j- 1] ^"j u\ K^{KR) 



+S3 



dv cosh I 



dv cosh 



-/? + ^ 



v] K,{KR) 



di/cosh ( P + iiO-j^] v) K^KR) 
Using the Kontorovich-Levedev Transforms J6|, we obtain 



(29) 



^ = !|(exp{-/fi?cosh(/3 + z[0-(j-l)^])} 
-Kexp|-ii:i?cosh('-/3 + i 0-{j + l)^ )\ 



exp < —KR cosh ( /5 + 
+ exp \ —KR cosh ( —j3 + i 



+53 



■{ 



exp \ —KRcosh 



+ exp \ —KR cosh ( (3 



{-'' 



+ i 



- {j + 1) 

")} 



} 



^3 



)} 



Introducing fi from En. lf23ll into this expression, yields 



(30) 



* 



(expj- 



6 



-i?cos 



0' - 1) 



ifci?sin 



'-«-"!]} 



■ exp <^ — it cos 



-0 + 1)3. 



ikRsiii 



" + ')il} 



+S 



exp -j — i?cos 



r TTC 

!si 



-0 + 1)5 



-ikRsin 

' TTC 

■exp ■! — -i?cos 
6 



o + "f]} 



- (i - 1) 



ikRsin 



+S3 



exp I ^Rcoa 



TTC . 



} 



■'3 



ikRsin 

+ exp|— Ecos [^- J3 



'ifci?sin 



J =0,1,..., 5 



}]) 



(31) 



From Ea. ()17|l it is easily seen that this wave function is fully symmetric under 
the interchange of particles. The function is invariant under the addition of 7r/3 
to 0, together with the addition of one unit to j, which is what should be done to 
move from one region in the (p, 0)plane to its next counterclockwise neighbour. 
Remember that for each region there is a specific order of the particles. 

We can also see that the real part of each of the exponential arguments is 
negative, except when — {j ^ l)'''/3 = ±7r/2, that is on the lines 9 = 9j, 
where its value is zero. Thus, when R is large, the wave function is negligible 
except near the lines 9 — 9j. Note that only the first 4 terms give a significant 
contribution in the asymptotic region. We can be conclude that the form of 
the wave function is that of products of bound state functions, associated with 
two particles, with oscillatory functions, which describe the location of the third 
particle with respect to the 2 bound ones. Evaluating, then, the wave function 
for large values of R, its asymptotic form can be written as 



^iR,9') 



^{-Kcose'} 



^{-ikRsine'} , <^ {ifc_Rsin9'}\ 



where tt/6 < 9' = 9—{j — 1) |- < 7r/2, j = 0, 1, ..., 5. This asymptotic expression 
consists of a wave representing a two particles bound state multiplied by an 
incoming wave, together with an outgoing wave multiplied by S. 
>From the expression in En. lf22|l the matrix S can be written as 



sin(f +i^)cos(| -i(i) 
cos(f +i/?)sin(| - i0) 
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sinf 


+ sin( 


-2if3) 










sinf 


+ 2 sin 


{i[3) cos (i/3) 










sinf 


— 2 sin 


(i/?)cos(i/?)- 




In terms of K, 
















m 


sin 


m = - 


TTC ik 

6K K 


—inck 






6(7r2cV36-fc 


') 


Therefore 


s^ 


1 


— 36 (k/nc) — i 


(24/\/3) (fc/^c) 






1 


- 36 {k/ircf + i 


(24/x/3) (fc/^c) 





sin ^ — sin (— 2i/3) 

(32) 

(33) 

(34) 

-I- 1 y^'i/y 6) yKj-Kc) 

This is, precisely, the scattering matrix 

^ ^ [-1 ~ » (6V3/7rc) fc] [3 + » (6V3/7rc) fc] 
[3 + i (6\/3/7rc) fc] [- 1 + i (6 VS/vrc) fc] ' 

given as Eq.(61) in Ref. 52|- It corresponds to the symmetric S matrix calcu- 
lated for the specific process 2+1. 

The matrix iSs, see again II22|I . has the following form as a function of k 

53 = ' + ^^'^^"^^' -1 + ,: (6V3/.c) fc. (36) 

^3 multiplies the shorter ranged part of the exact wave function, that goes 
to zero when R goes to oo. 

Additional insight can be gained, by following the reasoning of McGuire 
[Tfl] . The scattering of 3 asymptotically free particles, to 3 also asymptotically 
free particles, requires 3 (successive) collisions, and yields the part of the wave 
function associated with the calculation of the S-matrix. Intermediate stages, 
associated with fewer collisions, give rise to the shorter ranged part of the wave 
functions. A similar reasoning holds for the 2 + 1 processes. 

In conclusion, we have shown that this integration contour, and the choice of 
Bessel functions, have imparted the correct asymptotic behaviour. Furthermore, 
we can deduce from the asymptotic expression that the coefficient S represents 
the S-matrix. 



Relation to Adiabatic Theory 

The eigenfunctions of the following eigenvalue equation [l^ form a complete 
set of orthogonal hyperspherical adiabatic basis functions. Changing, a bit, the 
usual notation: 



+ ^)-;|c(e) + A.(i?') 



^ (^ + 7 ) - ^^ (^) + A. (i?') BM R') = 0, (37) 
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where R' , a real parameter in this equation, is held fixed; k — 0,6,12,... and 
Ak [R') -^ (k^ — l/4)/i?'^ as the interaction is turned off. The unnormalized 
eigenfunctions 

BJe-R')= cos 



where j is an integer such that 16* — j-3| < ^ and q^ satisfies 

7ri?'c 



(38) 



g^tan 



(g<Z«) = ^. (39) 



In the adiabatic approach the parameter R' is identified with the hyper radius 
R. For E <Q there is only one open channel, labeled by k = 0. For large R, the 
channel function is concentrated along the lines defined hy 9 = 6j . Accordingly it 
can represent the two-body bound state. For E > 0, there is an infinite number 
of open channels, labeled by the successive numbers k, equal and greater than 
zero. They describe, asymptotically, three free particles or a two-body bound 
state, together with a free particle. 

Note that the function p{i^), En. lfTsll . is a real function if, and only if, i' 
takes on values along the imaginary or the real axis, see Fig. 2. It can be 
seen that the pseudo-Sturmian function defined in Eq.Q coincides, apart from 
normalization constants, with the lowest adiabatic function Bq{0; R') when v = 
qo is an imaginary number and p{h') = R'{qQ). Also, li v = q^ are in the 
real intervals (3 + [k — 6] , 9 + [k — 6]) with k = 6, 12, ..., then p(j/) — R'{qK.) 
and the pseudo-Sturmian functions become equal, except for the normalization 
constants, to the adiabatic eigenfunctions Bi^{6;R'). 

Thus, in the case of the example considered in this paper, that is E < and 
the 2-1-1 system, the integral Ea. ljl2ll . along the imaginary axis in the complex 
v plane, can be written in terms of the lowest adiabatic function as 

^{R,e)= f diyA{iy)Bo (9; R'iiy)) Z^{KR), (40) 

where v runs from —ioo to ioo . The most important contribution of the adia- 
batic functions to the integral, at large R, comes from the lines 9 = 6j, where 
two of the particles are joined. When these adiabatic functions are multiplied 
by the appropriate Bessels functions, their linear combination (Eq. (I4fl|l ) should 
have the correct asymptotic behaviour, and will represent a two-body bound 
system in the colliding with a third particle. 



Conclusions and Outlook 

We have shown that the integral representation approach within the hyperspher- 
ical context, when applied to McGuire's model, offers a reliable tool to study the 
collisional dynamics of the 3-body system. We have obtained several interesting 
results, namely: 

-An exact solution to the corresponding Schrodinger equation. 
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Figure 2: Plot of the pseudo-Sturmian eigenvalue p{iy). In (a) we plot p{iy) as 
a function of v'^. In (b) the plot of (a) is rotated and flipped to give v^ as a 
function of p. For p positive, v'^ = A{p)p'^ + j. 



p(v) 



' I ' 7^ I •^^ • I ' I 
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-A closed form for the angular basis for this system, the pseudo-Sturmian 
functions. 

-A recurrence relation for the coefficients, in the expansion of the wave func- 
tion in terms of the free-particle basis. 

-The S-matrix, obtained directly from the solution of the recurrence relation. 

-The relation of the present approach to the traditional adiabatic approach. 

-The relation of the present solution to the known plane wave exact solution. 

The simpHcity of the approach as compared with the adiabatic one, promises 
to be very useful in extending it to more compHcated situations, Hke the system 
with different masses and systems with more particles, currently under research, 
or systems in three dimensions modeled by ZRP potentials. In the last case, the 
method can be applied to a wide kind of systems to obtain asymptotic solutions 
which can be matched to solutions obtained with methods like the R-matrix 
one, simplifying substantially the calculations. 
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APPENDIX A: Pseudo-Sturmian Functions 



Fixing u, the general solution for the eigenvalue equation 
92 
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vio) 






^{0), 



(41) 



with 9j — {2j + 1) 7r/6, can be written as the free angular wave solution if{9) = 
D^ cos [u {9 — 7j)] , j = 0, 1, ..., 5, provided that it be continuous through the 
boundary lines 9 — 9j , that is. 



cos [ly {9 J - 7j+i)] == cos [ly {9, - -j,)] , 
and satisfies the boundary conditions 



(42) 



lim 



e,+C 



s,-C 



d9 



d9^ ^ 

"^ (=0 



^(9) 



V{0) = 0, 



(43) 
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with j — 0, 1, ..., 5 . D^, which does not depend on j for the symmetric solution, 
determined by normaHzing the wave function [ig!. The requirement of continuity 
leads to the conditions 7^- +7j+i = 29j — {j + [j + 1]) n/S or to 7j_i = 7j. The 
second condition does not satisfy lll3|l so we shall use the first one, which can 
be written as 7j — j n/3. Now to focus on Eq.(E}. Continuity impHes that 
the integral of the second term gives zero. For each j, the first and the third 
terms give 



lim 



C^o va^(^°^'^[^-^^+i])e=fl,+. 



— (cosz^[6l-7j] 



-pC'^)-^ cosi/[6lj -7j] = 0. 



If we select a symmetric solution, then 



d d 

— cos {i^ [9j +i- 7j+i]) = — cos (i/ [9j - (^ - 7j]) 



(44) 



(45) 



and taking the limit in Eg . (J44II . we obtain the desired form of the boundary 
condition 



lim 



1 



d 



cos{iy[0-jj]) 



TTC 



(46) 



p{u)cos{v[9--fj])d9 
where j = 0, 1, ..., 5. Calculating the derivative and the limit in Eq. (lifill yields 



6 /^ / N 

— V tan j/TT/o = p {v) 



(47) 



We conclude that cos [j^ (0 — J7r/3)] , j — 0,1,. ..,5, satisfies Ea. (j41|l . pro- 
vided that p{i^) satisfies En. l|17jl . 

APPENDIX B: Derivation of the plane wave representation in 
terms of cartesian coordiantes 

For the 2 + 1 system, and aside from an ultimate normalization, the incoming 
wave function from En. ljMll can be written in terms of cartesian coordinates, as 



i^' 



exp 



TTC I Xi — X2 



6 V V2 
V6 



(j-1) 



(j-1) 



, xi +a;2 - 2x3 
—ik I — cos 



V6 



(j-1)' 
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Xi — X2 



V2 



■ sm 



-«-"i 



+ exp 



vrc / xi — X2 



V2 



O' + l) 



a;i + X2 — 2x3 . 
H ;= sm 

/Xi +X2 - 2X3 

— zfc = COS 



O' + i) 



V6 



(.? + !)■ 



Xl — X2 



V2 



0- + 1) o 



(48) 



Evaluating the trigonometric functions for j = in the above expression, the 
argument of the first exponential function takes the form 



2 , / . TTC . 1 , 

—;=k Xl + —I — ■= + i^-=k 
\/6 V 6\/2 V6 

1 



i^=k xa 



6\/2 ^/6 
By labeling the particle wave numbers as in p?j . 



2^2 



(49) 



fc2 = 



4fc, 



6\/2 ^/6 

TTC 



it. 



6V2 \/6 



(50) 



the incoming wave takes the form 

V"* = exp{-i(fc3Xi + fc2a;2 + A:iX3)}^^g 
+ exp {i (fc2Xi + fc3X2 + fcia;3)}j^o 

The outgoing wave for j — can be obtained from the incoming one by substi- 
tuting k by —k, which in turns means interchanging ki ^ ^2 and inverting the 
sign of the whole argument within all exponentials, that is. 



S exp {-i (fcixi + fc3X2 + k2X3)}j^Q 
+ exp {i (fcaxi + A;iX2 + A:2X3)}^^g 

The wave associated to the factor ^3 can be written as 

exp {-i (fcixi + ^2X2 + ^3X3)1^-^0 
+ exp {i (^2X1 + fciX2 + k3X3)}j^Q 
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The above results correspond to the sector j = in the (p, 9) plane, in which the 
order of particles is given by 0:2 < Xs < xi. The waves in different sectors can be 
obtained by the appropriate permutation of the set of coordinates {xi,X2,X3}. 
The completely symmetric wave plane may then be written as 



'j=]p 



V" = Sp I exp [-i {k^xi + k2X2 + kix^) 

+ exp [i {k2Xi + k3X2 + kix^)]-^^^ | 

+S{p) jcxp [-i {kixi + fcgxa + ^22:3)]^=^^ 

+ exp [i {ksxi + kiX2 + k2X3)]j^j^ | 

+53 (p) jcxp [-i {kixi + k2X2 + k3X3)]j^-^ 

+ exp [i {k2Xi + kiX2 + k^x^)].^-^ | (51) 

where the sum runs over all permutations of the set {xi,X2,X3}. 
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